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The ability to protect quantum information from the effect of noise is one of the major goals of 
quantum information processing. In this article, we study limitations on the asymptotic stability 
of quantum information stored in passive JV-qubit systems. We consider the effect of small imper- 
fections in the implementation of the protecting Hamiltonian in the form of perturbations or weak 
coupling to a ground state environment. We prove that, regardless of the protecting Hamiltonian, 
there exists a perturbed evolution that necessitates a final error correcting step when the state of the 
memory is read. Such an error correction step is shown to require a finite error threshold, the lack 
thereof being exemplified by the 3D compass model. We go on to present explicit weak Hamiltonian 
perturbations which destroy the logical information stored in the 2D toric code in a time 0(log(iV)). 

PACS numbers: 03.67.Lx, 03.67.Pp 
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I. INTRODUCTION 

Quantum information processing promises exciting 
new capabilities for a host of computational [T] [5] [3] and 
cryptographic [H[S] tasks, if only we can fabricate devices 
that take advantage of the subtle and very fragile effects 
of quantum mechanics. The theory of quantum error- 
correcting codes (QECCs) and fault-tolerance [5] [7J [5] H] 
assure that this fragility can be overcome at a logical level 
once an error rate per element below a certain threshold is 
achieved. However, providing a scalable physical imple- 
mentation of computational elements with the required 
degree of precision and control has proven to be a task 
of extreme difficulty Thus, one might hope to design 
superior fault-tolerant components whose robustness is 
enforced in a more natural way at a physical level. 

A first step in this daunting task is to concentrate 
not on universal quantum computation, but on one sub- 
protocol within this; the storage of quantum informa- 
tion. Thus, the aim is to find systems naturally assuring 
the stability of quantum information, just like magnetic 
domains in a hard disk provide stable storage of classi- 
cal information. The quest for such a passive quantum 
memory was pioneered by Kitaev [TU], who introduced 
the toric code as the first many body protecting Hamilto- 
nian. The promising conjunction of properties shown by 
his proposal has fueled a search, which is yet to provide 
a definitive result. 

For families of protecting Hamiltonians, such as Ki- 
taev's toric code [TO] [TTJ, a constant energy gap 7 sep- 
arates the degenerate ground space, used for encoding, 
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from low energy excited states. Furthermore, the stabi- 
lizer representation of these Hamiltonians naturally asso- 
ciates it with a QECC, which permits an error threshold 
without the use of concatenation [TT]. A perturbation 
theoretic expansion of local errors V in the Hamiltonian 
must then cancel to orders increasing with the distance 
of the associated QECC. Thus, within the region of va- 
lidity of degenerate perturbation theory, the degeneracy 
of the ground state space is only split by an amount that 
is exponentially small in the distance of the code. In 
turn, it takes an exponentially long time for the splitting 
to implement logical rotations on the perturbed ground 
space (e.g. a phase gate). 

However, such perturbation theoretic results must be 
applied with caution. The most important limitation 
probably arises from the fact that they deal with a closed 
quantum system whereas actual noise may be better 
modeled by perturbative coupling to an environment. 
Furthermore, the rigorous range of validity of pertur- 
bation theory [T^], \\V\\ < 7/2, is extremely restrictive 
when considering perturbations with an extensive opera- 
tor norm. Even if local observables can be adapted for to 
a high degree of accuracy, as shown by Hastings and Wen 
[IB] , the global eigenstates of the system may change and 
become very different. Within our understanding, the 
possibility of adapting encoding and decoding protocols 
relies on the perturbation being characterized 1 . However, 
most of our work will focus on the relevant scenario of an 
unknown perturbation, where unperturbed (unadapted) 
encoding and decoding protocols are used. 

Recently, Alicki et. al. have presented results support- 
ing the instability of quantum memories based on Ki- 
taev's 2D toric code [TS] and the stability of its 4D ver- 
sion [16] when coupled to a sufficiently cold thermal en- 
vironment. Chesi et al. |17j have made progress in pro- 
viding a general expression giving a lower bound for the 
lifetime of encoded information. The approach taken in 
these articles is thermodynamic in nature and has the ad- 
vantage of allowing the derivation of positive results. A 
weak coupling Markovian approximation to an environ- 
ment at thermal equilibrium is assumed, thus neglecting 
any memory effects from the environment. In a previ- 
ous article [T5] , we considered a Hamiltonian system sub- 
ject to independent depolarizing noise (corresponding to 
the high temperature limit of the above approach) and 
proved that O(logiV) is the optimal survival time for a 
logical qubit stored inside N physical qubits. 

Our current approach directly deals with Hamiltonian 
perturbations and environment couplings without going 
through a Markovian approximation for the environment. 
A comparative advantage of our approach is the capabil- 
ity of exactly dealing with certain weak but finite per- 
turbations and couplings, and providing restricted no go 



1 A possible exception to this is given by proposals of adiabatic 
state preparation 114) . 



results. 

We consider the effect of relatively weak yet unknown 
perturbations of an N qubit local protecting Hamiltonian 
and coupling to an ancillary environment starting out in 
its ground state. We show that as the number N of phys- 
ical subsystems used grows, it is impossible to immunize 
a quantum subspace against such noise by means of local 
protecting Hamiltonians only. We further show that if 
one wishes to recover the quantum state by means of an 
error correction procedure, the QECC used must have 
some finite error threshold in order to guarantee a high 
fidelity; this result is applied to the 3D compass model 
which is shown not to have such a threshold. In the case 
of the 2D toric code, we propose Hamiltonian perturba- 
tions capable of destroying encoded information after a 
time proportional to log(iV), suggesting that some form 
of macroscopic energy barrier may be necessary. Weak fi- 
nite range Hamiltonian perturbations are then presented 
which destroy classical information encoded into the 2D 
Ising model; in this case interactions involving a large, 
yet iV independent, number of qubits are required. Fi- 
nally, we consider time dependent Hamiltonian pertur- 
bations and coupling to an ancillary environment with 
a high energy density; here we provide constructions il- 
lustrating how these more powerful models may easily 
introduce logical errors in constant time into informa- 
tion protected by all stabilizer Hamiltonians, and certain 
generalizations. 



A. Noise model motivation 

A prerequisite to assess protecting Hamiltonians is a 
precise definition of the noise model they will be expected 
to counter. Our aim is to understand the protection life- 
time they provide to (quantum) information as well as 
to identify the properties a good protecting Hamiltonian 
should have. In order to be able to make such predictions, 
we will study noise models admitting a mathematically 
tractable description while striving to keep our choices 
physically motivated. 

To falsify claims of protection against any possible 
noise of a certain class (such as weak local perturbations 
to the Hamiltonian) , it suffices to consider an adversarial 
choice within such a class. In such a noise model, dif- 
ferent perturbations and environments are not assigned 
probabilities; a perturbation is simply considered possi- 
ble if it adheres to certain conditions. Finally, we assume 
that the perturbation has not been characterized. This 
allows us to derive no-go, or limitation, results from the 
exact analysis of adversarially engineered noise instances. 

The most elementary way in which the Hamiltonian 
evolution of a closed system can be altered is by including 
a small perturbation V to the Hamiltonian H. A simple 
physical interpretation for such a perturbation is to as- 
sociate V to imperfections in the implementation of the 
ideal protecting Hamiltonian H. Furthermore, Hamilto- 
nian perturbations extending beyond the system under 
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experimental control are modeled by a weak coupling 
between the system and an environment. We focus on 
families of protecting Hamiltonians satisfying certain lo- 
cality and boundedness conditions, and naturally extend 
similar restrictions on the perturbations and couplings 
considered. 

Let us first introduce some definitions. A family of pro- 
tecting Hamiltonians {ij/v} is parametrized by a number 
iV which grows with the number of physical subsystems 
participating in H^. A Hamiltonian H is called "fc-local" 
when it can be represented as a sum 

JV 

8=1 

with at most fc physical subsystems participating in each 
interaction term T\. The interaction strength of a phys- 
ical subsystem s in a fc-local Hamiltonian H is given by 
the sum ^\ ||Tj|| of operator norms over those interaction 
terms Ti in which the physical subsystem s participates. 
A family of fc-local Hamiltonians is called " J-bounded" if, 
for every Hamiltonian Hn in the family, the largest inter- 
action strength among the physical subsystems involved 
is no greater than J. Finally, a family of Hamiltonians 
will be D-dimensional if the physical subsystems involved 
can be arranged into a D-dimensional square lattice, such 
that all interaction terms are kept geometrically local. 

We will concentrate on families of fc-local, J-bounded 
protecting Hamiltonians, with J > 0, and fc, J ~ 0(1). 
Furthermore, the specific Hamiltonians treated in this ar- 
ticle admit an embedding into 2, 3 or 4 spatial dimensions 
and we may assume such embeddings also when dealing 
with generic protecting Hamiltonians. 

The families of Hamiltonian perturbations {Vjv} which 
we will consider will be J-bounded, with the strength 
J small in comparison to J. The perturbations will be 
taken to be fc-local, with fc possibly different, and even 
larger, than fc. This allows, for example, taking into con- 
sideration undesired higher order terms which may arise 
from perturbation theory gadgets |19j . Allowed pertur- 
bations should also admit a geometrically local interpre- 
tation under the same arrangement of subsystems as the 
protecting Hamiltonian. 

When considering coupling to an environment, an ad- 
ditional set of physical subsystems will be included as 
the environment state. A family of local environment 
Hamiltonians {H N } will be defined on these additional 
subsystems. The coupling between system and environ- 
ment will be given by a family of weak local Hamiltonian 
perturbations vjf E ^ , acting on both system and environ- 
ment. 

H N = H^®lW+I^®HW + V jf E) (2) 

Finally, it should be possible to incorporate the addi- 
tional physical subsystems from the environment while 
preserving the number of spatial dimensions required for 
the Hamiltonian. To simplify notation, the sub-index N 
shall in general be dropped. 



The engineering of fc-body interactions is increasingly 
difficult as fc grows [T9j [20]. This is why we limit our 
study to families of fc-local Hamiltonians (i.e. fc indepen- 
dent of N). It is under such criteria that we exclude 
proposals such as quantum concatenated-code Hamilto- 
nians [21 , for which the required degree of interactions 
would grow algebraically with the number of qubits. 

The J-bounded condition guarantees that the rate of 
change for local observables remain bounded. Further- 
more, this condition is strictly weaker than the rigor- 
ous requirements stated for the simulation through the 
use of perturbation theory gadgets [19]. There, constant 
bounds are imposed both on the norm of each interac- 
tion as well as on the number of interactions in which 
each subsystem participates. The J-bounded condition 
also leaves out systems with long range interactions, as, 
for those systems, the total interaction strength of indi- 
vidual physical subsystems diverges as the system size 
grows. Such long range interacting systems are physi- 
cally relevant, and may lead to protecting Hamiltonian 
proposals [2"2l |2"5] . However, we abstain from treating 
such models for which our notion of weak perturbation 
seems inappropriate. 

Each physical subsystem may independently be sub- 
ject to control imprecision. Such is the case for weak 
unaccounted "magnetic field" acting on every compo- 
nent of the system or a weak coupling of each component 
to an independent environment. Thus, relevant physical 
scenarios involve perturbations with extensive operator 
norm (i.e. scaling with the number of subsystems). The 
J-bounded condition encapsulates these scenarios and 
seems to better describe what we understand by a weak 
perturbation. 

Finally, it is expected that scalable physical implemen- 
tations should be mapped to at most three spatial dimen- 
sions. This would rule out the 4D toric code Hamiltonian 
, a proposal which was otherwise shown to provide in- 
creasing protection against weak local coupling to a suf- 
ficiently cold thermal bath [16] . As would occur with an 
actual physical embedding, we expect that the perturba- 
tions considered may be included into the same geomet- 
rical picture as the protecting Hamiltonian they affect. 



B. Outline of results 

In the following sections, we analyze the problem of 
obtaining increased protection for quantum information 
by means of an encoding and a protecting Hamiltonian 
acting on an increasing number of subsystems. We con- 
sider the effect of adversarial noise models consisting of 
local Hamiltonian perturbations and/or a weakly coupled 
environment. The aim is to examine the assumptions 
and limitations of memory schemes based on Hamiltonian 
protection with a growing number of physical subsystems 
as quantified by the survival time of stored information. 

We prove in complete generality that if the read-out 
process does not incorporate a recovery procedure then 
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information retrieval will be unreliable. The figure of 
merit considered here is S(t) = tr (|^>(0)) (ip(0) \pifj), the 
overlap between initial and evolved state after a constant 
time t. For arbitrary protecting Hamiltonians we provide 
a completely general construction involving a weakly cou- 
pled environment starting in its ground state (Sec. II B ) 
which yields an exponentially small (in N) upper bound 
on S(t) after a constant time. For gapped Hamiltonians, 
a proof proceeding without reference to an environment 
(Appendix |A| can provide an upper bound to the time 
averaged overlap which is close to \ . 

From this point on, protecting Hamiltonians are con- 
sidered together with a recovery operation 7Z, applied on 
read-out, thus providing a more robust figure of merit 
Sn(t) = tr (p(0)TZ(p(t))) . A similar weak coupling con- 
struction shows that information content of the 3D com- 
pass model [53] can be destroyed in constant time by a 
zero temperature environment (Sec. Ill), despite of its 
error correcting mechanism 1Z. This provides a direct, 
negative, answer to the open question of whether the 
3D compass model is a self-correcting quantum memory. 
From a broader perspective, the structure of our proof 
strongly suggests that the underlying QECC defining the 
recovery operation 1Z must have a strictly positive error 
threshold. 

We continue by considering the effect of Hamiltonian 
perturbations on the 2D toric code [10 . The recovery 
mechanism 1Z is then taken as the composition of a fixed 
syndrome measurement followed by a correction opera- 
tion pairing the detected anyons. It is shown (Sec. IV) 
that, although the underlying QECC has an error thresh- 
old, it is not protected against combinations of unknown 
weak local Hamiltonian perturbations, even after a fi- 
nal round of error correction is considered. Our claim is 
based on adversarial weak local perturbations that are 
capable of destroying the stored information in a time 
logarithmic in N. This is stronger than previous results 
[25) in that, the noise model requires no interaction with 
the environment and the information is destroyed expo- 
nentially faster. 

In a similar manner, we consider perturbations on 
the 2D Ising Hamiltonian (Sec. [Vf, which is often 
used as an example of self-correcting classical memory. 
Here, Hamiltonian perturbations may transform (classi- 
cal) code states into an ambiguous state in constant time. 
While the number k of bodies in perturbation terms is 
required to grow as the overall perturbation strength de- 
creases, it shows no dependence on the size N of the sys- 
tem. In this model, any sequence of local errors connect- 
ing the two classical code states must go through states 
with a macroscopic amount of extra energy, showing that 
this property alone is not sufficient to give protection. 

Beyond the Hamiltonian perturbation model and cou- 
pling to a ground state environment, we consider more 
aggressive noise models (Sec. VI) in which the environ- 
ment can introduce large amounts of energy. The mod- 
els considered are time dependent Hamiltonian perturba- 
tions and weak Hamiltonian coupling to an environment 



starting in a high energy state. For such noise models, 
even the information storage capabilities of the 4D toric 
code, a proposal shown to be thermodynamically stable, 
are completely destroyed. 

One might expect that the results presented here are 
not limited to the task of designing a quantum mem- 
ory. Rather, they tell us about the difficulty of keeping 
a state and its time evolution confined within a specific 
subspace of the system, under the effect of Hamiltonian 
noise. Such considerations arise in other settings, such 
as in the models of adiabatic and topological quantum 
computation. We will outline some of these connections 
in Sec. EH] 



II. NECESSITY OF ERROR CORRECTION ON 
READ-OUT 

We start our examination of passive quantum mem- 
ories by proving that a final step of error correction is 
necessary, regardless of the choice of encoding and pro- 
tecting Hamiltonian. What we mean by this, is that if 
such a recovery step is not included on read-out, there 
is no way of guaranteeing a high fidelity (close to 1) be- 
tween initial and final states for more than a constant 
amount of time. 

First, we propose weak local perturbations showing an 
exponentially decreasing overlap between perturbed and 
unperturbed eigenstates. For general local Hamiltonians 
and states, we consider a local coupling V of the system 
with a 7-bounded environment initialized in its ground 
state. Averaging over such couplings V, we are able to 
derive an exponentially small upper bound (S(tf))y < 
[1 — sin 2 (2e)] N for the overlap between initial and evolved 
system states at a time tf = —. 



A. Eigenstate susceptibility to perturbations 

Consider a fc-local Hamiltonian H , decomposable into 
fc-body interaction terms as described in (JT|). We 
choose a perturbation V such that the initial and final 
Hamiltonians arc related by a composition of local uni- 
tary transformations, 



H = H + V = UHU\ withW 



N 



1 = 1 



(3) 



where P; are normalized local Hermitian operators. Tak- 
ing this definition, V can be written as 



N 



(4) 



and thus is also fc-local. Furthermore, if 2fce 1, it is 
justified to call V a perturbation with respect to H, since 
all terms are small with respect to those of H. 
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Degeneracies in H are assumed to be infinitesimally 
lifted to ensure uniquely denned eigenvectors. The over- 
lap between eigenvectors \ipi) of H and the perturbed 
eigenvectors U\ipi) is then given by Fu — \ (ipi\U\ipi) | 2 . 

By averaging over all possible directions Pi, we effec- 
tively obtain an independent qubit depolarization. 

J U^i^UdPt ...dP N = AfgdViX^I) (5) 

Here, A\(p) = \p + (1 — A)| is the qubit depolarizing 
channel and A(e) — 1 — |sin 2 (e). We may then denote 
{F)u as average of the overlap Fy over all local rotations 
IA having a given strength e. This average is expressed 
in terms of the depolarizing channel as 

(F) u = (^\Af^m(^\)\^). (6) 

A result of King [2B], known as multiplicativity of the 
maximum output p-norm for the depolarizing channels, 
states that 

max||Ar(|0)(0|)|| p < (^max||A A (|0)(0|)|| p ^. (7) 

For qubit subsystems and for p — oo, Eq. ^ bounds 
the overlap of Af N (\<f>){<p\) with any single pure state, 
leading to 

(F)u < (^f = (i-i^Kf ' (8) 

Not only does this imply the existence of specific rota- 
tions such that F becomes exponentially small as the 
number of subsystems N grows, but that this is true for 
most rotations IA. While this is already known under the 
name of Anderson's orthogonality catastrophe (see, for 
example [27J HE]), we re-derive it for completeness and as 
an opportunity to introduce techniques needed through- 
out the rest of the paper. 



B. State evolution in coupled Hamiltonians 

In this section, we consider a weak Hamiltonian pertur- 
bation coupling the system to a "cold" environment. The 
environment is assumed to start in its ground state, cor- 
responding to a cold environment assumption. Averag- 
ing over a specific family of such perturbations instances 
V , an exponentially small bound on the overlap between 
the initial state and the evolved state is obtained. This 
bound, (S(tf)) v < [1 — sin 2 (2£)/3] Ar , is obtained after a 
constant evolution time tf = ^ , inversely proportional to 
the strength of the environment Hamiltonian. 

Suppose that we start with a state \ipo) "protected" by 
an N qubit system Hamiltonian Hg. We can introduce 
a simple environment, composed of 2N qubits, each of 



which starts in its ground state, |0), and which is defined 
by its Hamiltonian 

H = ® + g, H (E) (9) 

N 

h {e) = TEli+Xi+L^-looxoof). (10) 

i=l 

When necessary, we take the supraindices (S), (El) and 
(E2) to denote the system, the first, and second compo- 
nents of the environment respectively. While both en- 
vironment components will interact with the system, it 
is the presence of both which will allow a simple inter- 
pretation of the induced decoherence as a probabilistic 
application of local errors. 

We again use the trick of considering a perturbed 
Hamiltonian H — UHIA^ which results from the weak 
local rotations IA = ®, = i Uj of the decoupled Hamilto- 
nian H. The rotation elements will involve both system 

and environment components, Uj = e %e i j , where 

(s) 

the operators Pj are taken to be Pauli-like operators 
on site j of the system. 

The perturbation V = UHW — H must be decompos- 
able into small local terms. Such a decomposition for V 
is given in terms of the decomposition H^ s > = 
into at most fc-body terms. Each perturbation term 
Vi = UTiW — Ti has an operator norm no greater than 
2e/c||Ti|| and involves up to 2fc-body interactions 2 . The 
perturbation required to rotate the environment Hamil- 
tonian terms involve at most 3-body terms and a total 
norm bounded by 2e^. 

The initial state, |?/'o)|00) ,glJV will thus evolve into 
e~ mHU ^ \i) )\QG)® N . The survival probability is then 
S(t) = (ipo\ps{t)\vpo), wnere Ps(t) = tr B (p(t)), and 

p(t) =Ue- itH tf |Vo)(V>o | ® |00)(00|® Ar We itH W t . (11) 

Here, U may be explicitly decomposed as 

U = expiie^ Pj s) ® Xj si) ) 

= EpCos(£) JV - w (P)(isin(£)) l "(P)P^ ) ®X( E1) , 

(12) 

where p denotes a binary vector indicating the sites on 
which rotations are applied in Pp and w(p) is the 
weight of the bit string p (number of non identity factors 

in /' P S ). 

Now consider a time tf — - such that e~ lt,H trans- 
forms components of the environment from 1 10) to |11), 
while leaving components in state 1 00} unaltered. At such 



2 For ek <C 1, a further decomposition of such terms can be pro- 
vided in which subterms involving k + b bodies are of strength 
0(e b ), guaranteeing that the strength of terms decays exponen- 
tially with the number of bodies involved. 
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a time 1 1 , substituting U into expressions (111 allows ex- 
plicitly tracing over the environment to yield 



Ps 



(*/) = E cos2 ( £ ) 



2iV-tu(p)-ti;(q) s j n 2 ( £ ) u, (p)+ u, ( c i) 



P v e-^f H{S) P^){MP^ t ' H{S) P v 



(13) 

Thus, ps {if) may be considered as the density matrix re- 
sulting from the independent probabilistic application of 
the local unitary rotations prescribed by U on |V , o)(V , o|i 
followed by the evolution under the unperturbed system 
Hamiltonian, followed by a second round of random ap- 
plication of the local rotations prescribed by U. Defining 



£rAp) = P R P R + C 1 - p)p> ( u ) (S(t')) t 

and the Hamiltonian evolution 

W t (p)=e-^ CS V* (S \ (15) 

we can take p = sin 2 (e) and define p v %rt := 
(&>iLi£Pi;pj IV'oXV'ol, so that we can express S(tf) as 

S(t f ) = tr (p mrt e- lt ' Hs Pmr te ltfHs ) . (16) 

This is the overlap between a density matrix and its own 
unitary evolution, and can be upper bounded by 



S(t f ) < tr (p 2 mrt ) 



(17) 



In turn, using the fact that Pj are Pauli-like operators we 
may rewrite it as 



N 



S(t f ) < (Vol \^$£p i<2p{ i- P ) 1 (IVoXVoDlVo)- (18) 

Averaging over the Pauli-like operators, we obtain 

(S(t f )) v < (^\Af { N e) (\^)^o\Mo), (19) 

which is the overlap at time 1 1 averaged over the proposed 
family of weak perturbative couplings. Here A\(p) is 
again the depolarizing channel, and X(e) = is 1— §2p(l — 
p) . Using Eq. ^ , we obtain 



(S(t f )) v <[l--p(l-p)] N , 



(20) 



which by substituting p for sin 2 (e) yields 

(S(t f )) v <ll-sin 2 (2e)/3} N . (21) 

By averaging over different possible weak couplings, we 
obtain an overlap between initial and evolved states 
which is exponentially decreasing in N. 

The norm 7 of Hamiltonian terms in the environment 
should be bounded, since it is in part these terms which 
are rotated by U to introduce a weak coupling between 
system and environment. Thus, the proposed evolution 



time tf = — is constant. Furthermore, if one considers 
an environment of N semi-infinite chains of coupled two 
level systems (such as Heisenberg chains), it is possible 
to ensure that the overlap with the initial state is small 
for all times larger than t ~ -, rather than have the 
recurrences that arise from the discrete spectra of the 
described model. 

In appendix |A|, the evolution of an unperturbed 
eigenstate is considered under the effect of pure Hamil- 
tonian perturbations (no environment). In this case , a 
constant rate of change in the system state is guaranteed 
by an energy gap 7 in the system Hamiltonian. If an 
initial state belongs to an energy band separated from 
the rest of the Hilbert space by such an energy gap 7, 
we provide an upper bound on the time averaged overlap 



'£[0,4] 



< 



1 

2jt 



between initial and evolved state. 



Discussion 



We were able to show after a constant time tf, an expo- 
nentially large degradation of the overlap S(tf) in terms 
of N. If we are to find a benchmark by which to evaluate 
a memory scheme, we expect that the memory improves 
as more resources are dedicated to its implementation. It 
is now clear that many-body quantum states are inher- 
ently unstable with respect to the uncorrected overlap, 
i.e. S(t) is not the appropriate benchmark. 

The fact that we may count on N physical subsystems 
to implement a quantum memory should not exclude us- 
ing only one of them and ignoring whatever noisy evolu- 
tion is affecting the others. This corresponds to consid- 
ering an overlap reduced to the relevant subsystem and 
not on the whole state. Already such a simple idea guar- 
antees that information storage quality is non-decreasing 
with N. 

One may further generalize this by realizing that the 
relevant subsystem need not correspond to an actual 
physical subsystem. This corresponds to preparing code 
states of a QECC and performing a single round of error 
correction 1Z on read-out. The resulting benchmark is 
Sn(t) = tr (p(0)TZ(p(t))) , quantifying the quality of the 
information extracted from the evolved state and not of 
the state itself. Robust error corrected logical observ- 
ables 3 can analogously be defined via the extension pro- 
vided from the code space to the whole Hilbert space by 
the recovery operation 1Z. In the following sections, we 
shall consider protecting Hamiltonians associated with 
such error correcting codes and robust logical observ- 
ables. 



3 Alicki et al. |16| use the name dressed observables, whereas 
Chesi et al. |17| use the self-explanatory name of error corrected 
logical operators, which we will adopt. 
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III. LIMITATIONS OF THE 3D COMPASS 
MODEL 

A desirable property for a quantum memory is that 
any sequence of local operators mapping between dif- 
ferent logical code-states should have energy penalties 
which grow with the system size. It has been shown that 
this happens for schemes in four dimensions, such as the 
4D toric code [IT]. Seeking to provide such an exam- 
ple in three spatial dimensions, Bacon proposed the 3D 
compass model [24] . a scheme based on subsystem er- 
ror correcting codes [22] and requiring only 2-body near- 
est neighbor interactions. Furthermore, mean field argu- 
ments suggest that this model might show such an in- 
creasingly large energetic barrier. 

However, we will show that the zero temperature (lo- 
cal, but non-Markovian) environment construction of the 
previous section is capable of giving a false read-out from 
the code after constant time. We argue that this is due 
to the choice of the prescribed error correction protocol 
by showing that the same flaw is present for the 4D toric 
code if one uses a similar error correction technique. We 
show that this failure is a general feature of all quantum 
memory protocols for which the underlying quantum er- 
ror correcting code does not have a local error threshold 
(is unable to handle errors on a constant fraction of the 
sites) . 

For the 3D compass model, quantum information is 
first encoded into states of a 2-local Hamiltonian H$ de- 
fined on a TV x N x N arrangement of two level systems 
(where N is an odd number). 



N N-X 



Hs — — A E E (Xi : ijXi + i^j + Xij t jXij + i u 



i,j=i i=i 



(22) 



This is not a stabilizer code but a subsystem code. This 
means that a recovery operation need not correct certain 
errors which have no effect on the logical observables, and 
information may be preserved even if the recovered state 
is different. 

First, note that pairs of planes of operators Z{ = 
11,., z ij,l z i,j,l . i and Xi = H,. ; A'/ ,.,.V;. commute 
with the Hamiltonian H for all I. This also holds for 
logical operators, which consist of products along a sin- 
gle plane ZeZ, = ]J Z i:J j and X = Xi = X^ hJ 
operators respectively, for an arbitrarily chosen I. In the 
ground space, the choice of I is irrelevant, since the op- 
erators Z[,X[ all have +1 eigenvalues. Provided N is 
odd, X and Z anti-commute, giving a qubit algebra. In 
the presence of errors (outside of the ground space) , the 
error corrected logical observables will be defined as the 
majority vote among plane observables 



Z ec = maj,Z ; X ec = ma^Xi 



(23) 



where maj stands for a majority vote among the ±1 eigen- 
valued commuting operators. Measuring all pairs of ad- 
jacent planes Zi allows a majority vote error correction 



scheme to be performed on the value of the Zi plane ob- 
servables without extracting whether the corrected state 
yields +1 or —1 values for all such planes. 

Considering a perturbation on the system plus an en- 
vironment, as presented in section |il B By explicitly de- 
veloping the final expectation values for the observables 
of interest (tr [X ec ps{tf)^j and tr (X ec ps(tf)^J ), it can be 
seen that the information stored in the code will not be 
reliable after a time tf. For this, we can pick up from the 



evolved state of the system in Eq. ( 16 ) 



Ps(t f ) 




l^oXV-ol- (24) 



Since all plane observables (A; and Zi) of a given type 
mutually commute, and also do so with the Hamilto- 
nian, we can independently consider the probability of 
each plane observable having suffered a flip. If the Pi in 
Eq. ( 14 1 are taken to be single X or Z rotations, they will 



anticommute with overlaping Z\ or X\ plane observables 
respectively, changing their value upon an odd number 
of applications. Taking the Pi to be simple Z operators, 
the probability of flipping the value of an Xi plane ob- 



servables by applying &)f =1 £p;. P once is given by 



i<N 



Ppla 



^ cos^ v ~-^(£)sin 2j (e) 



^2N 2 -2i 
i£odd 

l-cos Ar2 (2e 



N 



(25) 



which is exponentially close to 1/2. Since all observ- 
ables involved commute with the system Hamiltonian, 
the probability of observing any result configuration will 
be preserved by H.t f . Finally, a second round of errors 

&)iLi£Pi,p wm again flip the observed value for each 
plane with a probability p p i anc *. The final independent 
probability of flipping the value of each plane is 

- 1 - cos2Af2 (2e) 

£*plane — ^Pplanc*l-L Ppl&nc* ) ~ ^ • \^) 

The proposed correction scheme is equivalent to a ma- 
jority voting among such planes. Thus, if more than half 
the planes suffer such an error, the majority vote will 
fail. The probability for incorrectly measuring the error 
corrected logical observable X ec on read-out is then 



E Ppla„c(l-Ppla„c) W -*(^). (27) 

Given that |[1 — cos 27V (2e)] < p p i a nc < \i we have that 

(28) 



Plog 



i[l-iYcos 2Ar2 (2e)] <p logic < \. 



Assuming e to be a small constant independent of N, 
the probability pi og icai will exponentially approach 1 /2 
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for large N. We conclude that the encoding is not robust 
against the error model posed by local coupling to a cold 
adversarial environment. 

The problem lies in the error correction mechanism 
rather than the protecting Hamiltonian itself. This be- 
comes apparent if one applies a similar analysis to the 
4D Toric code. There, the suggested noise model does 
not present a problem, since the usual error correction 
[TT] of the 4D toric code has an error threshold. That 
is, provided the probability of per-site error, p = sin 2 e is 
below this threshold, there exist error correction criteria 
which succeed with a probability approaching 1 exponen- 
tially with N . On the other hand, we could consider a 
majority voting version of error correction in this setting, 
where we measure hyperplanes of X operators, and apply 
a majority vote to choose the correct result. In this case, 
an analysis completely analogous to that of Bacon's 3D 
compass code would hold proving such a read-out tech- 
nique unreliable. We conclude that the error correction 
procedure of the compass code does not allow sufficient 
resolution to use any potentially topological properties of 
the encoded quantum information. 

The errors introduced by <2)j =1 £p;,p are sufficiently 
general to suggest a necessary criterion for Hamiltonian 
protection of information from weak coupling to a cold 
environment. Even if we consider the first round of er- 
rors and the Hamiltonian evolution as part of the encod- 
ing procedure, information should still be able to with- 
stand the probabilistic application of arbitrary local er- 
rors. This means that the information, either quantum 
or classical, should be encoded in such a way as to pro- 
vide a finite error threshold in the thermodynamic (large 
N) limit. Of course, in practice, it will be desirable to 
have a code with a fault-tolerance threshold such that 
when we try to implement the round of error correction, 
faulty operations can be compensated for. 



IV. LIMITATIONS OF THE 2D TORIC CODE 

We will now show how local Hamiltonian perturba- 
tions are capable of introducing uncorrectable errors in 
Kitaev's 2D toric code Hamiltonian. The introduction of 
such errors will strongly rely on the lack of string tension 
on the toric code, suggesting a macroscopic energy bar- 
rier may be a necessary requirement. A brief introduction 
to the toric code Hamiltonian is provided in appendix |B| 
and is recommended to the unfamiliar reader. 

Logical operations in the 2D toric code can be real- 
ized by creating a pair of anyons, propagating them so 
as to complete a non-trivial loop, and finally annihilating 
them. It is roughly such a scheme that will be followed 
by the perturbations we develop here. Repeating tech- 
niques from section[TTJ we may consider the initial state as 
containing a superposition of local errors which arc inter- 
preted as neighboring anyon pairs. A perturbation con- 
struction due to Kay |25j allows the deterministic prop- 
agation of such anyon pairs along predefined adversar- 



ial paths on the lattice. Syndrome measurement allows 
restricting to a probabilistic picture where error strings 
corresponding to anyon propagation paths are present 
with a predefined probability. It is finally the recovery 
procedure which may possibly complete these errors into 
logical operations by selecting an incorrect anyon match- 
ing. 

A family of weak local perturbations capable of prob- 
abilistically introducing distant anyon pairs will first be 
presented. As before, the initial state \i/j(0)) is assumed 
to be a ground state of the unperturbed Hamiltonian H, 
in this case an N x N toric code as in Eq. (Bl I. After 



a time tt proportional to the maximum desired anyon 
propagation distance D, unperturbed syndrome read-out 
on \4>(tf)) will probabilistically detect distant (as well as 
local) anyon pairs. Our construction will then be applied 
to produce a simple set of O(N) distance anyons such 
that no syndrome based error correction may be reliably 
applied. Later, shorter yet more elaborate anyon prop- 
agation paths will require explicit analysis of the error 
correcting probability of different anyon pairing proto- 
cols. In this context, we find weak Hamiltonian pertur- 
bations are capable of introducing logical errors with a 
large probability (~ ^) in a time tf logarithmic in the 
system size N. 



A. Probabilistic introduction of distant anyons 

Kay [25] showed that local errors (anyons) in the 2D 
toric code, and other local stabilizer Hamiltonians lack- 
ing string tension, can be propagated into logical errors 
corresponding to almost complete loop operators by a lo- 
cal Hamiltonian perturbation P. While in his work the 
initial presence of the anyons was assumed, here, anyons 
will be introduced with a certain amplitude by a gener- 
alization of the Hamiltonian perturbation P. 

Consider introducing perturbations of the form V = 
U{H+P)U^ —H , where U = ® 4 Ui decomposes into weak 
local unitary rotations, and P is, as in |25j . a weak lo- 
cal perturbation capable of deterministically propagating 
anyons in a given time tf. The perturbed Hamiltonian 



H = U(H + P)U i 



(29) 



induces a time evolution which can be written as 



= e-^I^O)) = Ue- u ( H+p W\il>(0)). (30) 
In this context, U and P are chosen such that: 

1. Neighboring vertex anyon pairs are created by U\ 
with a certain small amplitude O(e), by applying 
weak Z rotations on connecting edges. 

2. Each of the anyons is deterministically propagated 
by P along a predefined path. Thus, local exci- 
tation pairs become strings of errors defining new 
positions for the anyon pair. 
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3. Finally, U is unable to remove both the anyons cre- 
ated by W after at least one of them has been prop- 
agated. Moreover, if none were present, U creates 
an anyon pair with amplitude O(e). 

Propagation paths for each anyon are not allowed to 
overlap but are otherwise completely independent. The 
propagation of the i-th anyon along its path £W may be 
attributed to a specific component Pj of P = P. L . In 
turn, each component Pj admits a decomposition 



Pi 



j'=i 



(i) 



T, 



(31) 



(i) 



in terms of local interaction terms T p ^ q , where 
the anyon locations along the path £. As in [5S] 



are 
the 

scalar coefficients J- are chosen to implement a perfect 
state transfer [301 EU EH] and each term T Pjg implement 
a swap among vertex anyons on p and g. If p and g are 
neighboring vertices, T pq is defined as 



T Pj g = Z s frac{l - A p A q )2, 



(32) 



where A p and A q are the vertex stabilizer operators cor- 
responding to p and q respectively (appendix B]) and Z s 
is a Z rotation on physical site s corresponding to the 
edge connecting p and q. Furthermore, by allowing p 
and q to be next nearest neighbors, it is possible to have 
crossing anyon paths £(*' without having them over- 
lap in the anyon locations used. If vertices p and q are 
not neighbors, the same effect is obtained by substituting 
Z s in (32 1 for a tensor product of Z operators along an 



edge path pq from p to q, 



Tn 



Z, 



(1 - A p A q ) 



(33) 



s£ pq 



The distance D is the maximum number of steps 
among the different anyon propagation paths D = 
maxi \£^\. It will be taken as D = N/2 — 1 in sec- 



tion (IVB) and as D = O(logTV) in section (IVDl. Fix 



■(i) 



ing the strength of perturbation terms in P as J- 

~ \J | + 1 — j) allows the perturbation P to remain 
e-bounded while allowing simultaneous perfect anyon 
transfer in a time tf = ^r. Similarly to previous sec- 
tions, by taking the rotations Uj — e leZj as e weak, 
the final perturbation required V will also be composed 
of O(e) strength interactions involving at most 8 bodies 
each. 

The quantum state before measurement at time tt is 



given in Eq. ( 30 1 . Expanding W from Uj = e leZj , we get 



\i/>(t f )) = Ue- i ( p+H )*t (g)(coselj- - isineZj)\ip(0)), 

(34) 

where the index j ranges over sites of non trivial action 
for U. The state | -0(0)) is a ground space eigenstate of 



H , and assuming the locations j on which U acts are non 
neighboring, each Zj will increase the energy respect to 
H by 7. Furthermore since the energy of a state respects 
to H depends only on anyon number and P is anyon 
number preserving, we have [H, P] = 0, allowing us to 
write 

|V>(*/)> =Ue- iPt f(Q(cosElj -ie-^f aneZj)]^)). 

(35) 

Since all the propagations in P commute and correspond 
to exact transfer of each anyon created by W precisely 
at time tj , we may write 



\tp(t f )) = WjI(cosel 



ie irytj sine 



W(0)) 



(36) 



where ££^' is the path given by the union of {j} and the 
two propagation paths £w + > and of P correspond- 
ing to the each of the two anyons created by Zj. By 
expanding U, we obtain 



mf))=U 



cos 2 el-ie ntf sine cose 6^) Z; 



+i cos e sin eZ, + sin 2 ee" i7t/ Z, 6d Z t 



(37) 



|V(0)>. 



The state \ip(tf)) described by Eq. (37) corresponds to 



a coherent quantum superposition of applying different 
error paths. For such unitary evolutions, initially orthog- 
onal states will remain orthogonal and thus fully distin- 
guishable. However, there are at least two mechanisms 
which lead us to consider a mixed density matrix as the 
final state. The first, is due to the fact that the actual 
perturbation applied is not known, and can for instance 
be taken probabilistically among the family of perturba- 
tions described. The second, is unperturbed syndrome 
measurement M, which is the first step of a quantum 
error correction procedure to recover the initial state. 

Syndrome measurement M will probabilistically 
project the state \ip(tf)} into a subspace consistent with 
a fixed anyon distribution. This is the first step of the 
recovery operation 1Z — CoJ\4, the sequential application 
of unperturbed syndrome measurement M. followed by a 
syndrome dependent correction operation C. Analysis of 
different correction strategies C need only focus on the 
resulting mixed state A4(\ip(tf))(tp(tf)\). Since for any 
anyon configuration there is at most one combination of 
operators yielding it in Eq. (37 1, the state \ip(tf)}(ip(tf)\ 
is reduced to a probabilistic application of these operators 
on |^(0))(V>(0)|. Again, taking £ R , p {p) = pRpR+(l-p)p, 
one may verify that 



M(m I ))m f )\) = o j £ Zj , P £< i 



Zi , P (Mtf))m f )\) 

(38) 
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Case Probability 

A (1-p) 2 

B p(l-p) 

C (l-p) P 

D v 2 



Errors 
...III ...III .. 
. . ZZZ . . . ZIZ . 
...IZI ...III .. 
. . ZIZ . . . ZIZ . 



Anyons 
. . . 0000 . . . 0000 . 
...0000... 0110. 
...0110... 0000. 
...0110... 0110. 



FIG. 1: Two hollow dots indicate positions where a pair of 
vertex anyons may be created by ~l0 and/or by U with prob- 
ability p. Anyons created by U' are propagated by P along 
the darkening path. A table is provided indicating the proba- 
bility of possible error configurations and their corresponding 
syndrome observables (1 (0) representing anyon presence (ab- 
sence)). 



with p — sin 2 e. Note that the order of application is arbi- 
trary, since the £r, p superoperators commute. Thus, one 
may consider independent probabilities p for observing 
each anyon pair created by W and propagated by P (or 
unpropagated anyon pairs created by U). Hence, when 
instantiating the Hamiltonian perturbation described on 
a certain set of anyon propagation paths, one need only 
deal with the independent probabilities of measuring 
propagated and unpropagated anyon pairs. 



B. Simple error loops in O(N) time 

The aim of this subsection is to provide a simple en- 
semble of perturbations, employing the above construc- 
tion in such a way that resulting anyon configurations 
are provably ambiguous, by which we mean that a sin- 
gle anyon configuration could have, with equal likelihood, 
originated from logically inequivalent errors. This means 
that for such configurations, the anyon pairing recovery 
procedure C can do no better than guessing, and will 
complete a logical error with a 50% probability for any 
possible choice of C. 

Let us first consider weakly perturbing only in the 
vicinity of a single row. The joint effect of many such per- 
turbations, will then be shown to produce further degra- 
dation of stored information. So, W introduces a single 
Z error (neighboring anyon pair) on a physical start site 
s of the row with probability p. The paths for the per- 
turbation P are chosen such that both of the produced 
anyons propagate along the row in opposite directions 
up to final neighboring locations which are diametrically 
opposite s (see figure [TJ. For s weak perturbations, this 
requires no more than 0(N/e) time. Finally, with prob- 
ability p, U may introduce an error at site s or counter 
an clement of the propagated error chain. As can be seen 
from the figure, if the anyon introduction site s is chosen 
uniformly at random, there are observable anyon config- 
urations which occur with probability 2p(l — p), which 
are completely ambiguous (e.g. cases B and C are indis- 
tinguishable under exchange of initial site s). However, 
if such a syndrome is measured, the correction proto- 



col has a 50% chance of completing a horizontal Z loop 
on the lattice, which is equivalent to applying a com- 
pletely dcphasing channel on one of the encoded qubits 
with probability 2p — 2p 2 . 

By applying such a perturbation family to i < N rows 
of the lattice, the probability of not having such a logi- 
cally dephasing action take place becomes (1 — 2p+2p 2 ) 1 , 
which may be made arbitrarily small for large N (i.e. 
an odd number of horizontal Z loops is completed with 
a probability exponentially close to 1/2). Completely 
analogous string like perturbations exist for any of four 
logical operators defining the 2-qubit algebra associated 
to the ground space. Again, by simultaneously consider- 
ing such perturbations on a sufficiently large set of par- 
allel lines these operators too will be completed with a 
probability exponentially close to 1/2. Furthermore, by 
allowing anyons to hop directly to next nearest neighbors 
(i.e. Eq. ( 33 1 ) , it becomes possible to simultaneously in- 
troduce perpendicular yet commuting loop operations as 
a result of anyon removal. 

Simultaneously introducing the four logical operators 
independently with probability exponentially close to |, 
would yield a state exponentially close to a maximal mix- 
ture over the code space. Our proof requires terms from 
different perturbation paths to commute, indicating a 
possible obstacle to achieving this. In practice however, 
given that different anyons follow roughly ballistic trajec- 
tories with a relatively small spread, this does not pose 
an issue. In appendix [C] we show how it is possible to se- 
lect the set of anyon trajectories in the perturbation such 
that the order of anyon crossing is well defined (exponen- 
tially well in N). In turn, this implies an exponentially 
small deviation from the result of performing such anyon 
propagations in order, resulting in a state exponentially 
close to a maximal mixture on the four dimensional code- 
space. 



C. Localization in 2D stabilizer codes 

In the perturbations constructed to introduce logical 
errors in the toric code, there is a strong use of the en- 
ergy degeneracy of subspaces with the same number of 
anyons. The strengths of the different stabilizer terms in 
the 2D toric code manifest as strengths of local magnetic 
fields in the effective Hamiltonian of the propagation [25j . 
However, having exactly the same strength for all local 
Hamiltonian terms is not an essential feature of the 2D 
toric code or of stabilizer Hamiltonians in general. 

In the unperturbed picture of stabilizer Hamiltonians, 
excitations are completely localized. However, when dif- 
ferent excitations live in a degenerate energy space, per- 
turbations may be very effective at propagating them. 
With the hope of obtaining some localization with re- 
spect to anyon propagation terms, different stabilizer 
term strengths may be randomly chosen from some range 

Tupper ^ Tlowcr ^ 0. 

However, a Hamiltonian perturbation may "smooth" 



11 



this distribution to take on, at random, only a finite num- 
ber of discrete energy values, separated by e, the strength 
of the perturbation. The number of such possible values 
is given by [ 7 " pp ° r ~ 7l ° wer ] , which is therefore also the av- 
erage spacing between sites of the same energy. Hence, 
by selecting propagation terms of a similar size, the hop- 
ping scheme can route around any defects, and give a 
logical error. Thus, the argument is unable to guaran- 
tee protection against any constant sized perturbation. 
Nevertheless, it may be that the perturbation terms nec- 
essary to break the code should involve a larger number 
of bodies, which would definitely be an improvement. 

In the case of the 2D toric code , and all other 2D 
local stabilizer Hamiltonians [33l [34], there are always 
logical operations with string like support. This means 
that, albeit with some possible energetic smoothening, 



the scheme presented in section IV B can be adapted to 
introduce logical errors in arbitrary 2D stabilizer Hamil- 
tonians, meaning that the asymptotic lifetime which 2D 
N x N stabilizer codes may guarantee against weak local 
perturbations cannot be more than O(N). 



D. Logical errors in 0(log(iV)) time 

In the previous section, we gave a rigorous upper 
bound of O(N) on the information lifetime of the toric 
code. This bound coincides with the one provided by 
Kay in |25j . which required initial anyons in the system 
to be introduced by an unspecified environment. In this 
subsection, we provide an exponentially tighter bound 
by concentrating on specific choices for error correction 
protocols. We argue that it is possible for a Hamilto- 
nian perturbation to introduce ambiguous distributions 
of anyon configurations in a time logarithmic in TV, i.e. af- 
ter a time tf ~ 0(log N), error correction succeeds with 
probability not much higher than 1/2. Figure [2] schemat- 
ically presents one such perturbation, indicating where 
anyon pairs should be introduced, and paths Pk along 
which they should propagate. The fact that the trajecto- 
ries have only simple crossings allows them to be imple- 
mented by weak local Hamiltonian perturbation terms 
involving at most 8 bodies, as obtained from Eq. (33 1, 



with p, q being next nearest neighbors. Furthermore, the 
trajectory length is no more than twice the distance at 
which the anyon pair is finally separated. 

The length of anyon propagation trajectories is 85 + 4, 
where 



S = 



In TV 

~2p 



(39) 



and each has 2S simple crossings with other trajecto- 
ries. The time required to perform such a propagation 
by fixed strength local perturbations is proportional to 
S (i.e. logarithmic in N). 

A relevant property of such a perturbation is that 
anyons observed when performing unperturbed error cor- 
rection after an evolution time tf are always collinear. 




FIG. 2: Anyon pairs corresponding to each thick red edge 
may be created by . After a time t f , the right anyon from 
each pair introduced will be propagated a distance 4S + 2 to 
the right introducing Z errors along the darkening paths. Fi- 
nally, U acting on the same red segment may move an unprop- 
agated anyon one position to the right or create a neighboring 
anyon pair on it. The number of big steps (or equivalently 
of crossings) during the upward propagation is given by S, 
which in the case of the figure is 2. 



The anyon type and line direction may be chosen to co- 
incide with any of the logical operations, translating to 
the fact that any logical error may be introduced. This 
also has the desirable effect of simplifying the analysis of 
anyon matching criteria. There are only two logically in- 
equivalent anyon matchings on the line, which are the two 
perfect matchings in which each anyon is paired with one 
of its two nearest neighbors (i.e. right or left). The point 
is that one matching will be logically equivalent to the 
actual trajectories performed by the anyons, canceling 
any errors introduced, whereas the other will complete 
the actual paths into a logical error. A simple criterion 
to determine which case we are dealing with is to count 
how many times the actual trajectories, together with 
the anyon matching, cross a vertical line or any homo- 
logically equivalent curve. An odd number of crossings 
means that a logical error has been completed, whereas 
an even number of crossings means that the proposed 
pairing has been successful at error correcting. 

We study the success probability of two apparently rea- 
sonable matching criteria. The first minimizes the fur- 
thest distance among paired anyons. The second, for 
which a polynomial algorithm is known |35j . consists of 
minimizing the sum of distances among paired anyons. 
Proofs and numerics will be provided for the large N 
regime given by N ^> 45 + 2 which convey a high logical 
error rate. 



Anyon matching that minimizes Loo 

Let us first consider minimizing the furthest distance 
among paired anyons. This is the L^ norm of the vec- 
tor with components given by the individual distances 
among anyons paired by the matching. We will prove 
that the probability of introducing a particular logical 
error is close to 1/2 by considering two disjoint scenar- 
ios. The first is the very unlikely scenario in which, on 
syndrome measurement, two consecutive anyons are mea- 
sured at a distance > D (by consecutive, we mean no 
additional anyons were measured in the interval between 
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them). The second is composed of anyon distributions 
consistent with the measurement of a fixed pair of consec- 
utive anyons at a distance < D. For such distributions, 
the number of activated anyon paths passing completely 
over the fixed pair is shown to be odd with probability 
very close to 1/2. 

Let us first bound the probability of observing two con- 
secutive anyons at a distance greater than D in the syn- 
drome measurement for the evolved state. Given a fixed 
region of length D, at least [Z? / 4J different potential 
anyon paths start and end in it. Furthermore, assum- 
ing D < AS, the probabilities for not measuring anyons 
in this region are independent and are 1 — p for each end 
of an anyon path and (1—p) 2 for each start of an anyon 
path, since both U and ift could have created anyons in 
this case. The anyon-free region can begin in any of ./V lo- 
cations of the full loop. Thus, regardless of correlations, 
the probability of having D consecutive anyon-free sites 
is upper bounded by JV(1 — p) 3 L- D / 4 J , 

Assume now that a pair of consecutive anyons is mea- 
sured at a distance no greater than D. There are at least 
[(4S—D)/A\ potential anyon paths going over this region, 
each with independent probability p of being observed. 
On syndrome read-out, the number of such paths that is 
activated is odd with a probability approaching 1/2 at 
least as fast as |(1 ± (1 - 2p)L( 4S -- D )/ 4 J). Since the L m 
norm correction completes a logical error if the most dis- 
tant consecutive anyon pair is covered by an odd number 
of activated anyon paths, then by inserting S — \ 2p 1 
and D = \ 8 ] , we get a probability lower bound for 
logical errors which approaches 1/2 as 1/2(1 — iV -1 / 5 ). 



Anyon matching that minimizes L\ 

Let us now consider the anyon pairing criterion that 
minimizes the total sum of distances among paired 
anyons. Since all anyons are found on a loop of length N, 
this criterion will always choose a pairing with total dis- 
tance no greater than N/2. Thus it will successfully error 
correct if and only if the total distance of regions of the 
loop covered an odd number of times by observed anyons 
is no greater than N/2. By taking S to be l"^^! , we ex- 
pect to find roughly half of the sites flipped. To see this, 
note that, on average, each site is covered approximately 
Se times. Moreover, the probability of each site being 
covered an odd number of times is |[1 — (1 — 2p) s \. For 
small p and the chosen value of S, the average number of 
sites covered an odd number of times is approximated by 
T — JN' Furthermore, we expect the actual number of 
such sites to approximately follow a normal distribution 
around this value, which would imply that logical errors 
are completed with a probability of close to \. However, 
since the flipping of different nearby sites are highly cor- 
related events, it is not clear how to go about proving 
this. Instead, computer simulations (Fig. [3]) provide very 
strong numerical evidence. 
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FIG. 3: The average probability of error for L\ correction af- 
ter the system evolves for a time tf under the described Hamil- 
tonian perturbation. Here, anyon pairs arise, and evolve to 
distances of [20 ln(iV)] , with a probability of 10%, all collinear 
on a line of length N. Each point represents an average over 
10 6 random samples, with error bars representing the magni- 
tude of estimated statistical errors. 



E. Discussion 

We have proven that Hamiltonian perturbations can 
completely destroy the information stored in the 2D toric 
code in a time proportional to N . The only assump- 
tions are that the precise Hamiltonian perturbation is 
unknown, and that recovery begins by performing un- 
perturbed syndrome measurements. A simple family of 
Hamiltonian perturbations with associated probabilities, 
was used to justify that the introduction of logical er- 
rors in O(N) time is fully independent from the correc- 
tion protocol used. This approach remains applicable for 
arbitrary 2D stabilizer codes, even when the stabilizer 
terms are of uneven strength. 

Furthermore, we have argued that logical errors may 
be introduced by weak local perturbations in a time log- 
arithmic with the system size. In particular, two appar- 
ently reasonable anyon pairing schemes were shown to 
provide an unreliable recovery mechanism against weak 
local perturbations acting for 0(log N) time. 

A fully general proof, including all possible error cor- 
rection strategies based on syndrome measurement is cur- 
rently lacking for the 0(logN) error introduction. The 
generality of the O(N) construction is obtained by con- 
sidering a family of different perturbations which could 
produce the same syndrome outcomes through topologi- 
cally inequivalcnt error paths. It may be fruitful to apply 
such an approach for a general proof of logical errors pro- 
duced in O(log-ZV) time. 



V. LIMITATIONS OF THE 2D ISING MODEL 

In present day classical computers, magnetic domains 
are widely used to provide passive safekeeping of classi- 
cal information. The Ising model is usually used to elu- 
cidate the origin of such long lived magnetized states as 
a collective effect arising from microscopic local 2-body 
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interactions 

H Iaiag = -J^2z i Z j . (40) 

(id) 

In two and higher spatial dimensions, the nearest neigh- 
bor Ising model presents a finite temperature phase tran- 
sition between a disordered phase and an ordered mag- 
netized phase. However, it has long been known that 
such a system looses its asymptotic bistability under the 
bias produced by even the weakest of magnetic fields 

|35l El ESI 

(ij) J 

Such studies consider the dynamics of minority droplets 
in a 2D Ising model as given by phcnomcnological equa- 
tions or the Metropolis algorithm. For a Metropolis al- 
gorithm in which such a systematic magnetic field e is 
included, there is only one stable phase parallel to the 
field. The anti-parallel phase becomes metastable, with 
a lifetime exponential in J/e. Dependence of the infor- 
mation lifetime on lattice size N appears only for the 
small N < J/e and will thus not appear if one first takes 
the limit for large N. 

In this section we consider storing one bit of classi- 
cal information subject to a quantum evolution of a per- 
turbed 2D Ising Hamiltonian. The observable on which 
the classical bit is encoded is assumed to be the overall 
direction for magnetization Z = maj^Zj. The perturba- 
tions may conceptually be split into two parts, Z paral- 
lel magnetic fields which introduce additional degenera- 
cies to the Hamiltonian, and transverse magnetic fields 
or many body terms which couple the new ground states, 
introducing a hopping between them. The perturbation 
terms considered will not show support or intensity grow- 
ing with N, and they will be capable of introducing log- 
ical errors to the unperturbed logical observable Z in a 
time also independent of N. 

A. Hamiltonian perturbation proposal 

Consider dividing the N x N 2D periodic lattice with 
a chessboard pattern of squares of M x M spins where 
M > 4J/e max and J is, again, the nearest neighbor 
Ising coupling constant and e max is the greatest local 
perturbation strength one expects the Hamiltonian to 
protect information against. For simplicity, we assume 
N = 2nM, where n is an integer. Consider alternately 
introducing iieZj magnetic fields in the lattice site j be- 
longing to white/black squares of the chessboard pattern 
respectively. The value e is chosen homogeneously for 
each square such that the energy difference, 2eM 2 , from 
fully field parallel and anti-parallel configurations of each 
square exactly matches the maximum energy difference 
for border Ising terms 8MJ. For N 3> J/e such a per- 
turbation is always possible. 



The point is that now, the ground space of the sys- 
tem acquires a much higher degeneracy, i.e. between 
2 2n +1 — 1 and 2 4 ™ . Each black square could be fully 
magnetized parallel to its preferred field direction or par- 
allel to the direction of its four neighboring squares if it 
is opposite. By taking one spin variable for each square, 
the ground states may be identified with those of a 2D 
nxn anti-ferromagnetic Ising model with magnetic field. 
Three important ground states are the two fully magne- 
tized lattice configurations and the checkerboard config- 
uration in which all spins are fully aligned to their local 
magnetic field. The gap of these ground states with re- 
spect to low lying exited states is 2s, which is the energy 
penalty of flipping a corner lattice site of a square that 
is fully oriented in the direction of the field but anti- 
parallel to the two neighboring squares adjacent to the 
stated corner. 

A flipping term for each square of the chessboard 
should be of the form a^Xj, where the j is taken 
over all the M 2 sites in the square. Such terms can 
be introduced either on all black squares or all white 
squares. This would respectively couple one of the two 
fully magnetized configurations with the checkerboard 
configuration, achieving a full swap of state in a time 
<Aip = -. This evolution is exact when such M 2 -body 
terms of norm a are allowed, which implies that a proof 
of Hamiltonian stability will not only require assuming 
sufficiently weak perturbation terms but also a specific 
bound for the number of bodies on which such terms act. 

Let us now focus on the magnitude of a. This is the 
coefficient for a many body term, in which the size of 
the support scales like M 2 = (4 J/e) 2 , independent of 
N. One may consider obtaining such a term from the 
M 2 -th order degenerate perturbation theory expansion 
of fields of the form ei J^j Xj . For perturbation theory 
to be strictly valid, one needs M 2 e2 < e. Even then, this 
small magnitude must be taken to the M 2 -th power to 
obtain the first non vanishing expansion term. The time 
required to flip all spins in a plaquette is then propor- 
tional to: 




This expression has no dependence on N and the same 
perturbation can be introduced in all squares of a given 
color to yield a fixed flip time. Furthermore, we note that 
the state of those chessboard squares which are not per- 
turbed is fixed and may be traced out exactly. Hence, the 
second set of perturbations applied are fully independent 
and degenerate perturbation theory may be rigorously 
applied. 



B. Discussion 

Although the flip time shows no dependence on N , it 
grows faster than exponentially in terms of AJ/e. It may 
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well be that for magnetic domains, describable by such 
a 2D Ising Hamiltonian as Eq. (40 1, the ratio 4J/e is 



sufficiently large to provide a lifetime longer than would 
be experimentally verifiable. More importatntly, we are 
dealing with an extremely simplified model, with the par- 
ticularity of neglecting any long range interactions of ac- 
tual physical systems. 

The fact that the perturbation is unknown means that 
if such a checkerboard state is observed on read-out, in- 
formation is not recoverable. Such schemes may clearly 
be generalized to higher dimensions and to deforma- 
tions of the checkerboard pattern. The existence of such 
perturbations elucidates important limitations for state- 
ments one may formally prove about the classical mem- 
ory reliability of the Ising model, and therefore what con- 
clusions one might draw about the presence of a macro- 
scopic energy barrier (string tension) which the 2D Ising 
model certainly possesses. However, it is not clear that 
these arguments can be applied, for instance, to the 4D 
toric code since it is a feature of classical memories, but 
not quantum ones, that local fields can split degeneracies. 



A. Time- varying Perturbations 

When considering Hamiltonian perturbations, we as- 
sumed that we were unable to determine the new ground 
space due to the perturbation, and thus encoded in the 
original ground state space. One might consider an in- 
termediate setting where encoding can be achieved in the 
perturbed code-space, perhaps due to an adiabatic evolu- 
tion such as proposed by |14j , or by a precise determina- 
tion of, and compensation for, any perturbations present 
at the start of the storage time. However, in real experi- 
ments, stray fields etc. responsible for perturbations may 
fluctuate in time. Again, if one can track the changes in 
perturbations, the proof of Hastings and Wen |13j con- 
tinues to hold because Lieb-Robinson bounds apply to 
time-varying local Hamiltonians, and we can therefore 
adapt the final error correction step as well. Instead, we 
proceed assuming it is impossible to precisely learn this 
time variation. 



Adiabatically varying perturbations 



VI. AGGRESSIVE NOISE MODELS 

In previous sections, we explored the effects of Hamil- 
tonian perturbations on quantum memories, and partic- 
ularly on the 2D toric code. We also considered exam- 
ples of local cold environments perturbatively coupled to 
a system, as illustrated by Sees. |IIB| and |III| The only 
energy available in these scenarios was due to local per- 
turbations on the system plus environment. Intuitively, a 
small but constant energy density proportional to e was 
allowed. While this energy is potentially 0(N d ) for a d 
spatial dimension lattice of N d qubits, it is difficult to 
concentrate it in specific regions in order to generate log- 
ical errors. In comparison, stabilizer codes only require 
0(N d ~ 1 ) energy to implement a logical gate through lo- 
cal rotations. 

More aggressive noise models may locally introduce 
large amounts of energy into the system while keeping 
perturbation magnitudes weak. Such an example is pro- 
vided by weak yet time dependent Hamiltonian pertur- 
bations. These are relevant when one considers effective 
protecting Hamiltonians in the interaction picture [39] . 
Another possibility is to consider the weak coupling of 
the system to an environment which starts in a high en- 
ergy state. Noise constructions for these models shall be 
presented in this section. 

In calling such noise models aggressive, we convey the 
fact that we do not expect "reasonable" Hamiltonian pro- 
tection schemes to guarantee a long lifetime against such 
models. Thus, their study may help identify required re- 
strictions on the noise model in order to allow for prov- 
ably robust Hamiltonian protected memory models. Fur- 
thermore, it may provide insight regarding potentially 
fruitful proof techniques. 



One extreme case to consider is that the perturbation 
varies adiabatically, so that the system remains in its 
ground state space. If we do not apply error correc- 
tion, then we are concerned with how long it takes before 
the initial and final ground states have a small overlap. 
We shall assume that the original Hamiltonian H of N 
qubits has an energy gap 7, and we will consider the 
time- varying perturbation 



V = U{t)HW{t) - H 



where, as before, 



JV 



WW = 11' 



-iteXj/T 



and T is the total time of the evolution, i.e. small local ro- 
tations are gradually introduced. At any time < t < T, 
the effective Hamiltonian U(t)HW(t) has the same en- 
ergy gap as H, which means that the adiabatic condition 
is satisfied for T ~ 1 /poly (7). From previous considera- 



tions, Eqn. (A9l, we know that the overlap of the initial 



state IV'(O)) an d the evolved state, the ground state of 
the adiabatically perturbed Hamiltonian, have an aver- 
age overlap of no more than tr (Pq) (1 — | sm 2 (te/T)) N . 
For large N and small e, this means that the final over- 
lap is of the order tr (P ) exp(— 3e 4 N ) if a phase of error 
correction is not involved. 

When error correction is introduced to this scenario, 
this maps into the situation where our quantum mem- 
ory is initially encoded in the perturbed subspace, but 
decoding is using the original, unperturbed, error correc- 
tion strategy. In the specific instance of the perturbation 
Li(T), we find that X rotations are applied probabilis- 
tically on each site, and hence our QECC must have a 
superior error threshold. 
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Hastings and Wen [T3] reveal a similar interpretation 
holds for all possible perturbations (assuming the gap 
remains open as the perturbation is introduced) since all 
local terms are converted into quasi-local rotations. 

2. Rapidly oscillating perturbations 

Another extreme scenario is when perturbations are 
allowed to oscillate with arbitrary frequencies. A sim- 
ple construction shows that for stabilizer Hamiltonians, 
this allows the introduction of arbitrary errors in con- 
stant time. We expect that optimal control theory may 
provide the tools to generalize such results to arbitrary 
Hamiltonians. 

Consider a stabilizer Hamiltonian H and a logical er- 
ror to implement L = PmPm-i ■ ■ ■ P2P1, which is decom- 
posed into Pauli operators Pi on different sites. We then 
consider the time dependent Hamiltonian perturbation 

V(t) = e(t) e- lHot Pe lHot . (43) 

i 

This perturbation is weak if e(t) is sufficiently small. Fur- 
thermore, given that H is a stabilizer Hamiltonian, V(t) 
may be written as a sum of local terms (at least as local as 
the stabilizer operators). Finally, the dependence of e(t) 
on time, is to allow for e(0) = which makes the initial 
encoding equivalent for both perturbed and unperturbed 
Hamiltonians. 

The point of such a perturbation, is that it is possible 
to explicitly calculate the evolution of the system state 
in the interaction picture. 

|Vi(t)> = n ie -^/° £ (*') dt '|<Mo)) (44) 

This means that after a constant time tt such that 
I = Jq* e(t')dt' , the target operation L is perfectly im- 
plemented in the interaction picture. If |^(0)) is an eigen- 
state, then L is also implemented in the Schrodinger pic- 
ture, modulo a global phase. 

Taking L to be a logical operator of the stabilizer code 
used, this means that time dependent perturbations of 
sufficiently high frequency can destroy stored informa- 
tion in constant time. Here, sufficiently high frequency 
refers to having perturbation terms which oscillate with 
frequencies at least as high as those corresponding to lo- 
calized excitations. 



B. Stabilizer Hamiltonians and energetic 
environment 

In what follows, we consider a model in which an envi- 
ronment starts out in an arbitrarily energetic state. How- 
ever, the couplings between system and environment are 
required to remain small and local. 

For simplicity, we assume that the system is defined by 
a stabilizer Hamiltonian H$ and that it starts out in an 



eigenstate |V>o) of all stabilizer operators. We will con- 
sider a sequence of M Pauli operators on different sites 
L = PmPm-i ■ ■ P2P1 compounding to a logical oper- 
ation. In the case of translationally invariant stabilizer 
codes, explicit constructions for these operators are given 
in [34 . Finally, we may assume a code state \ipo), such 
that (^oWo) - 0. 

Motivated by the realization that, in order to intro- 
duce logical errors, we need to transfer some energy from 
the environment to the system, we choose a specific en- 
vironment Hamiltonian He = —Hs ( a t this point, the 
complex conjugate is unnecessary, but will become use- 
ful later). This means that all steps up in energy in the 
system correspond to an identical step down in energy 
in the environment. We start the environment state in 

m- 

In this scenario, the coupling 

M 

H SE = sY, P s,i®P*Ea (45) 

is enough to produce the logical error L in constant time 
To see this, consider the two states Ps,i\iMPe ilV'o) 
and Ps,iPs,i\ipo) Pe iPe^o) ■ Here, the subindex i is an 
arbitrary binary vector indicating which values of j a 
product of Psj (respectively P E • ) should be taken over. 
First of all, since we have assumed H$ is a stabilizer, and 
the Ps,i are Pauli operators, the aforementioned states 
are zero eigenstates of Hg £g> 1^ + I5 ® He- Further- 
more, the effective Hamiltonian for the perturbation term 
sPs,i <& P% i acting on the pair of states Ps,\\i>a)PE ilV'o) 
and P S ,iPsMo)P E ,i P EMl) is .i ust a matrix 

£ (? 0)' 

independent of i. This means that we can consider the 
action of the different Pg.i <8> P% i terms independently: 

e-' lHsEt \%) = ^e- iePs ><® p »A |Vo). 

Due to the effective Hamiltonian, P$, i | ipo ) P E i|*0o) ^ s 

mapped to Ps^PsA^o) p ea p e.M*o) in a time 7r /( 2e )- 
Thus, the effect of the entire perturbation is to rotate, in 
a time 7r/(2e) from |^ ) to L\ipo). 

Of course, this approach requires the environment 
state to have a very high initial energy, namely to start 
in one of its highest energy states. A refinement of this 
argument allows us to only change the sign of stabilizers 
in He which share support with L. For a local stabilizer 
code in d spatial dimensions consisting of N d qubits, it 
was shown [51] that there are logical operators L with 
support on k cx A^ -1 sites. The initial state of the en- 
vironment \iPe) — I "00 ) i s still an eigenvector of these 
stabilizers, with the same eigenvalues. Thus, coupling 
to an environment with an energy proportional to N 11 " 1 , 
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may also introduce logical errors in the same time. This 
means that the energy required from the environment 
per system qubit tends to as (compare to perturba- 
tions, which introduce an energy e per site). The catch 
however, is that the distribution of the energy in the ini- 
tial environment is highly specific and is in general very 
different from distributions that may be provided by low 
temperature thermal states. 

We conclude that no stabilizer Hamiltonian will be ca- 
pable of providing a guarantee for the logical integrity of 
stored information under the presence of an adversarial, 
weakly coupled, local environment. Further statistical 
assumptions such as energy distribution associated to a 
low temperature environment state need to be included 
in addition to the weak local coupling assumptions. 



C. Non-stabilizer Hamiltonians 

Stabilizer Hamiltonians are not the only possible can- 
didates for providing information protection, although 
they are particularly attractive because local errors re- 
main as local errors (not propagating or multiplying in 
the absence of perturbations). Let us now consider the 
more general case of distance preserving Hamiltonians 
i.e. ones which might not leave local errors perfectly lo- 
calized, but do not increase the distance of the error as 
defined by an error correcting code 4 . Using the same 
construction as in the previous section, we will show that 
weak coupling to an environment can also introduce the 
relevant logical error into a distance preserving Hamilto- 
nian (i.e. a logical operation converting between the most 
distant code states) for classical memories, by which we 
mean that one set of local errors becomes irrelevant, say 
Z errors, and the presence of a logical error on the clas- 
sical bit depends only on the local X errors present. The 
distance preserving assumption means that the number 
of X errors is preserved, [Hs,Y^i^i\ = 0. The maxi- 
mum distance between any 2 states is for the eigenstates 
1 0) i and 0^ |l)i, suggesting we should use these states 
for encoding. 

Similarly to the previous subsection, we introduce an 
environment, and a perturbative coupling between sys- 
tem and environment, 

H(e) = H s <8> 1 E + Is <8> H E +e^I, s (8 X^ E . (46) 

i 

The system Hamiltonian is weakly coupled to a "mir- 
ror" system He = —Hg. This perturbative coupling 



is responsible for the evolution of a mirrored state 
|0)® JV |0) lglAr , eigenstate of the unperturbed Hamiltonian 
H(0). 

To analyze the evolution, let us consider the action of 
the operators in terms of the eigenstates of Hs. Due 
to the commutation relation, there must be (^) eigen- 
states {IVVnj)} of Hg which are simultaneous eigenvec- 
tors of J2i %i with eigenvalue 2m — N. We can thus 
express the eigenvectors \ipm,j) of Hs in terms of the 
canonical basis as 

MW>= E «(m)u*i|0> 8JV = E ^m) td \i), 



«(i)= 



,(i)= 



(47) 

where i are binary vectors with m non zero components 
and |i) are the respective states from the canonical basis. 
The matrix a(m) is unitary as it relates two orthonormal 
bases of the same subspace. Define 

1 



E 



M 



(m) i> fc >W 



E a ( m )rj a ( TO )i,fcivvj)iv4,fc) 



E ivw>i$ 



From this, one obtains that 



(Hs 



l s ® H* s ) \m) = 0, 



(48) 



(49) 



implying that any non trivial evolution of |m) arises ex- 
clusively from the perturbative coupling and is given by 



H(e)\m) — eJ„ 



1) +eJ m+ i|m+l) 



(50) 



with J m = \J m{N + 1 — m). These are precisely the 
coefficients performing perfect state transfer between 
|0) = |0)® 2iv and \N) = |1)® 2JV in a constant time t = § 

isniisa. 

These results exclude the possibility of proving robust- 
ness against weak adversarial coupling to an arbitrarily 
initialized environment, even of many classical memories 
using the repetition code (such as Ising models). We 
learn that if the environment can provide enough energy, 
then even weak local couplings may be sufficient to pro- 
duce logical operations. This also motivates the desire 
to encode in the ground state space of the Hamiltonian 
since, were we to encode in a higher energy subspace, the 
environment needs less energy to cause destructive ef- 
fects. Alternatively, the mechanism presented here could 
present a useful way to implement gates on a memory 



4 As an aside, note that Bacon's 3D compass code 1241 is an exam- 
ple of a code where the errors do not remain in fixed positions, 
but preserve the values of the observables (since the observables 
commute with the Hamiltonian), which reminds us (Sec. |IIl| l that 
error correction is already likely to become much more problem- 
atic for these codes. 



VII. FURTHER APPLICATIONS 

Constructed perturbations and results presented in 
this article have focused on elucidating limitations of pas- 
sive quantum memories. However, our results may be 
recast in the following other scenarios. 



17 



Adiabatic Quantum Computation- The standard 
approach to adiabatic quantum computation consists of 
implementing an adiabatic evolution 

H(t) = f(t)Hi + (l-f(t))H f , (51) 

between Hamiltonians Hi and Hf, where /(0) = 1 and 
f(T) = 0. While the ground state of the initial Hamilto- 
nian Hi is expected to be readily prepared, the ground 
state of the final Hamiltonian Hj encodes the result of 
the desired quantum computation. An energy gap no less 
than 7 between ground state and excited states of H(t) 
is required for the duration of the adiabatic evolution. 

In this context, it is possible that Hamiltonian per- 
turbations could change the initial or final ground state, 
and maybe even close the gap during the Hamiltonian 
trajectory. For example, a time dependent perturbation 

V(t) = {UH{t)U^ - H(t)) , (52) 

with U defined as in Eq. (§, can make the perturbed ini- 
tial and final ground states almost orthogonal to the un- 
perturbed versions (see Eq. (|8|), while keeping the same 
gap as H(t). Even assuming the perturbed initial ground 
state is exactly prepared, only if the final state belongs to 
a code space with an error threshold, will it be possible 
to reliably recover the desired result, as in |VI A| 

Connections between adiabatic quantum computation 
and passive quantum memories can be expected to con- 
tinue into the regime where error correction is incorpo- 
rated, and future studies may better elucidate the issues 
involved in developing a fault-tolerant theory of adiabatic 
quantum computation |40j . 

Topological Quantum Computation- Difficulties 
in implementing quantum memories can also be related 
to some of the difficulties in implementing a topological 
quantum computation. In particular, in section |IV| we 
illustrated how constant Hamiltonian perturbations can 
create and propagate anyons in the 2D toric code. In 
the context of topological quantum computation, where 
gates are implemented through the braiding of anyons, 
the existence of perturbations capable of creating and 
propagating anyon pairs is at least equally disturbing as 
in the memory scenario. 

Quantum Simulations- One of the most interesting 
uses of a quantum computer is likely to be the simu- 
lation of other quantum systems. While one could ex- 
press these simulations in terms of the circuit model of 
quantum computation, and from there create a circuit- 
based theory of fault-tolerance for quantum simulation, 
it would be advantageous to understand how this could 
be implemented more directly, via the simulation of an 
encoded Hamiltonian. 

A logical first step would be to encode the state of 
each subsystem to be simulated into a quantum mem- 
ory. Thus, establishing when quantum memories exist, 
or when they fail, has implications in this case. One of 
the most commonly applied techniques in Hamiltonian 
simulation is that of the Trotter-Suzuki decomposition, 



where pulses of non-commuting Hamiltonians are com- 
bined into one effective Hamiltonian to some accuracy 
5. This inaccuracy may be treated as a time dependent 
Hamiltonian perturbation. Given the power such pertur- 
bations were shown to have, it is with great care that one 
should consider the use of passive quantum memories as 
elements for such quantum simulators. 



VIII. CONCLUSIONS 

In this paper, we have studied several constraints on 
the extent to which a many body Hamiltonian can be ex- 
pected to protect quantum information. First of all, we 
showed that gapped local Hamiltonians have cigenstates 
which are asymptotically unstable under local Hamilto- 
nian perturbations. This result, commonly referred to as 
Anderson's orthogonality catastrophe [27] shows that a 
gap is not sufficient to guarantee protection against errors 
[231 |H] . We proved that a weakly coupled cold environ- 
ments, can in constant time, alter the evolution of any 
quantum state leading to an exponentially small overlap 
between initial and final states. Taking these results to- 
gether, we conclude that quantum memory schemes must 
incorporate at least one final round of error correction. 

Similar perturbations show why the 3D compass model 
[24] , a self-correcting quantum memory proposal, is not 
reliable, revealing that the code and error correcting pro- 
cess must posses an error correcting threshold. Similar 
conclusions may be drawn in scenarios where informa- 
tion encoding and evolution is according to a perturbed 
Hamiltonian but read-out and decoding are not. 

Further explicit counterexamples illustrate that while 
we may, in some cases, be able to adapt for known per- 
turbations, arbitrary unknown perturbations can destroy 
the storage properties of codes such as the 2D Toric code 
in a time 0(log N). In this case, the proposed adversarial 
Hamiltonian perturbation heavily relies on the absence of 
a macroscopic energy barrier (it is possible to transform 
orthogonal encoded states via a sequence of local oper- 
ations while keeping intermediate states in a low energy 
subspace). By considering the 2D Ising model, we have 
argued that, in and of itself, a macroscopic energy barrier 
is not sufficient to protect against perturbations. 

Finally, we have considered strong noise models such 
as time varying Hamiltonian perturbations and weak 
coupling to an arbitrarily initialized environment. We 
showed that these noise models could apply logical trans- 
formations on information protected by stabilizer Hamil- 
tonians or distance preserving classical memories in con- 
stant time. We expect these result to provide insight into 
how one may prove properties of passive quantum mem- 
ories and under which assumptions. For instance, since 
such time- varying Hamiltonian perturbations can destroy 
the 4D toric code, then when trying to prove robustness 
against static perturbations, Licb-Robinson bounds are 
unlikely to be beneficial. 

Having proven a variety of limitations for quantum 
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memory models and elucidated some required conditions, 
the next step is to incorporate this deeper understand- 
ing into new designs for quantum memories. One ma- 
jor route is to establish a set of necessary and sufficient 
conditions under which a quantum memory is protected 
against unknown weak static perturbations. Under such 
a model, we may once again raise the question of whether 
good protecting Hamiltonians in two or three spatial di- 
mensions exist. Furthermore, one would hope to find sim- 
ilar conditions under an extended perturbation model al- 
lowing a perturbatively coupled local environment. Here, 
a central problem is to determine which physically real- 
istic assumption may be made on the environment such 
that positive results are still attainable (i.e. conditions 
on the initial state of the environment, such as it being 
prepared in its ground state). Finally, one may study 
the possibility of engineering an out of equilibrium en- 
vironment to provide additional protection to quantum 
information. 
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APPENDIX A: STATE EVOLUTION IN 
PERTURBED GAPPED HAMILTONIANS 



Energy gaps are considered as a positive feature for a 
protecting Hamiltonian, since they are expected to pro- 
vide an energetic barrier which an error process is re- 
quired to overcome. However, it will be shown that for 
sufficiently large N, the fidelity of the unperturbed eigen- 
states acquires an upper bound close to 1/2 after being 
evolved under the effect of a perturbed Hamiltonian for 
a time inversely proportional to the gap energy 7. 

If a system is perturbed, but we do not know the na- 
ture of the perturbation, the best strategy is, arguably, to 
continue using the unperturbed encoding (i.e. the eigen- 
states of the unperturbed Hamiltonian). The survival 
probability for an unperturbed eigenstate |?/>o) of H after 
evolution under a perturbed Hamiltonian H for a given 
time t (Eq. ([3])) is, without error correction 

S(t) = \(Me- U "m\ 2 = K^olWe-^lVo)! 2 , (Al) 

i.e. we can express S(t) as the overlap of W\ipo) with it- 
self under the evolution of the unperturbed Hamiltonian 
H. Furthermore, in terms of the eigenstate decomposi- 
tion 

Wt|^ )=5^a J #i) where Htyj) = E^j), (A2) 




FIG. 4: There is an energy gap 7 separating the eigenenergies 
corresponding to an exponentially small subspace Po from the 
energies of the Hamiltonian eigenstates giving rise to the rest 
of the Hilbert space. 



S(t) may be expanded as 

S(t) = J2 I 1 2 la/cosK^ - E )t\. (A3) 

Assume the initial state \4>o) belongs to an energy sub- 
space Po of H (i.e. (ipo\Po\ipo) = 1), and that H imposes 
an energetic gap 7 between the subspace P an d its or- 
thogonal subspace (see figure |4| . This allows the sum in 
Eq. (|A3l to be split as 



S(t) 



E 

\1>i),\*l>j)ePo 

E 

l*>,h/>i>0fb 

2 E 



M 2 K| 2 cos((^ - Ej)t) 



cti\ 2 \aj\ 2 cos((Ei - Ej)t) 



a,\ 2 |a J | 2 cos((£: i - Ej)t) 



(A4) 



hMefb,|V>3>g.Po 

We then define R, the UPoW subspace overlap of \ipo) as 

r = {^\up u^\iI> )= J2 H 2 - ( A5 ) 

Taking the time average (S(t'))t'E[o,t] — \ Jo S(t')dt', 
and noting that 



\Ek -Ej\> 7 => 



cos((£i - Ej)t')d£ 



< -, (A6) 



Poincare recurrences are averaged out, providing a bound 



(S(i')>t'e[o,t] < R 2 + (1 - Rf + ~R(l - R). (A7) 
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Although the bound in Eq. (A7) is minimized for R 



1/2, this does not imply that the smallest values for 
(S(t')) t i,z[ t ] are actually obtained for R = 1/2. 

A sufficient condition for the existence of a weak per- 
turbation yielding R = \ may now be obtained by means 
of continuity arguments. First, note that R depends con- 
tinuously on the parameter e appearing in the definition 
of the rotation U. and R = 1 for e = 0. This means that 
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if, for some eo > 0, we find that R < 1/2, then R must be 
equal to 1/2 for some smaller positive value < e < Sq. 

As in the previous subsection, we may take (R)u as 
an average of the overlap R over different directions of 
the rotation IA. An expression for {R)u, in terms of the 
depolarizing channel is given by 



(R)i 



MPoAf^d^-X^I) 



(A8) 



Including the dimension of the subspace Pq, the same 
bound as in Eq. (|| may be used, leading to 



(R)u <tr(P ) ( l-^sin 2 (e) 



N 



(A9) 



If the asymptotic growth of t r (Po) is slower than 

(l — |sin 2 (e)) , the bound (A9l will be exponentially 
decreasing with N. This means that for sufficiently large 
iV", and for most directions of rotation, there is some small 
rotation parameter e yielding R — 1/2. For the impor- 
tant case of small e and a constant dimension tr(Po), 
large N refers to N ~ 0(e~ 2 ). 

For those U leading to R = |, the time averaged sur- 
vival probability (<!5(i'))i'e[o,t] for the corresponding per- 
turbation may be bounded as 



{S{t')) v 



e[o,t] 



< 



1 



1 

2yt' 



(A10) 



We thus obtain that the overlap of initial encoded states 
and uncorrected evolved states will drop to values not 
much larger than | in a time inversely proportional to 
the gap 7. 



APPENDIX B: THE TORIC CODE 

Kitaev introduced the toric code [10] with the inten- 
tion of achieving reliable storage of quantum information 
at the physical level, as in classical stable storage, rather 
than by periodically performing explicit error correction 
procedures. He proposed that the Hamiltonian of the 
physical system being used to store the quantum infor- 
mation could, by its nature, make the information stable. 
His proposal consisted of a 2D system with non trivial 
topology (such as the surface of a torus) with a stabilizer 
Hamiltonian composed of local terms. Qubits could then 
be stored in the ground subspace with a degeneracy of 
4 9 , with g being the genus of the surface on which the 
physical qubits are located. 

In the toric code Hamiltonian, the physical qubits are 
located on the edges of a planar grid covering the 2D 
surface. For concrctcncss and simplicity, we shall restrict 
to the case were the surface is a torus and the grid is 
an N x N square lattice (i.e. 2N 2 physical qubits). The 
Hamiltonian is composed of commuting terms which are 
products of Pauli operators on different sites (it is a sta- 
bilizer Hamiltonian) . For each vertex s of the grid, there 
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FIG. 5: Each edge in the grid represents a physical qubit 
and opposite sides of the grid are identified by toric periodic 
boundary conditions. Typical plaquette and vertex operators 
are depicted near the center. Two vertical loop operators, X\ 
and i?2, which allow breaking the degeneracy are also pre- 
sented. One can take these to be the X and Z operators 
for the first and second logically encoded qubits respectively. 
The complementary (anticommuting) operators are given by 
analogous horizontal loops. 



is a star (or vertex) term A s = IljGstar(s) which is the 
product of X operators over all the qubits of edges reach- 
ing s. Analogously, for each face p of the grid, there is 
a plaquette (or face) term B p = EliebouadaryCp) Z ] which 
is the product of Z operators over all the qubits of edges 
surrounding the face p. Since each vertex and face have 



either or 2 common edges, the terms A s and B p always 
commute. Hence all terms of the toric code Hamiltonian 



H 



(Bl) 



commute, and may be simultaneously diagonalized. 
Since H B A S =I and ]J p B. p = I, there are only 2N 2 - 2 
independent binary quantum numbers asociated to these 
terms (stabilizer operators) and each valid configuration 
determines a subspace of dimension 4. Due to this, vi- 
olations of plaquette (vertex) conditions A s \rp) = 
(Bplip) — \tp}) always come in respective pairs. Following 
usual nomenclature, virtual particles called vertex (pla- 
quette) anyons are respectively associated to these exci- 
tations. The set of stabilizers may be completed with a 
pair of logical observables consisting of the product of Z 
(X) operators along non contractiblc loops on the lattice 
(dual lattice), which may not be expressed as a product 
of plaquette (star) terms as illustrated in figure [5] To- 
gether with the set of Hamiltonian stabilizers any com- 
muting pair of these four logical operators (Xi, X 2 , Z\ 
and Z2) uniquely determine the state. 

A stated prerequisite for using the toric code as a pro- 
tecting Hamiltonian is that the energy splitting of the 
ground space due to Hamiltonian perturbations should 
be small. This is argued through the use of degenerate 
perturbation theory and the fact that it only gives non- 
zero splitting when the order taken is at least the lat- 
tice width/height, claiming an exponential suppression 
of perturbations in the ground space. 
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The interaction terms in this Hamiltonian may be used 
as the syndrome measurements of an error correcting 
code with the desirable property that they are all geo- 
metrically local. Such codes provide a way of obtaining 
a fault tolerance threshold without requiring the use of 
concatenated quantum error correction. In this case, in- 
creasing the lattice size allows periodic measurements to 
suppress the effect of errors up to any desired accuracy 
[TT] provided the accumulated error probability between 
measurements is below a certain threshold. 

We will briefly review how error syndromes are inter- 
preted and corrected, making the simplifying assumption 
that the error syndromes are measured perfectly. These 
syndromes, i.e. measurements of the stabilizers, reveal 
the presence of any anyons on the lattice, but do not 
distinguish between them, so it is up to us to determine 
how these anyons should be paired up in order to annihi- 
late them. For each of the two kinds of anyon, the error 
correcting procedure will pair up the anyons and anni- 
hilate them by applying a connecting string of operators 
on them. If the connections performed and the actual 
origin of the anyons form topologically trivial loops (con- 
tractible loops), the error correction will have been suc- 
cessful. If however, the actual error pathways, together 
with the connections performed by the error correction 
procedure complete one, or an odd number of, non-trivial 
loops, then a logical error will have been implemented. 

Different criteria for pairing anyons may lead to log- 
ically different results. This is illustrated in figure [6] 
where two different criteria are used to pair up six anyons. 
In particular, if one of the criteria compensates the ac- 
tual error path, allowing recovery of the initial state, the 
other will complete the error path into an undesired log- 
ical operation. 

There are two correction protocols which we will con- 
sider, as they are expected to perform adequately when 
correcting a small proportion of randomly located errors. 
The first, which we refer to as Zi-EC, consists of min- 
imizing the sum of distances among paired anyons, for 
which there is a polynomial time algorithm [35 . The 
second, Zoo-EC minimizes the furthest distance among 
paired anyons. 



APPENDIX C: FULL DEPOLARIZATION OF 
THE TORIC CODE'S PROTECTED SUBSPACE 



In the main body of the paper (Sec. IV B), we gave a 
construction for a single logical error Xi , X 2 , Zi or Z 2 to 
be applied with a probability exponentially close to 50%, 
independent of the model used for error correction. This 
is not sufficient to show that we get full depolarization 
of the two-qubit subspace because it is not automatically 
clear that all 4 logical errors can be introduced simultane- 
ously in the same model; the problem being that crossing 
paths for anti-commuting operations do not necessarily 
have a well-defined phase, and the perfect state trans- 
fer operations can fail. Indeed, if two non commuting 



FIG. 6: Illustration of a possible configuration of three vertex 
anyon pairs (small circles). Segments indicate possible qubits 
where Z rotations could be introduced in order to remove 
the anyons. Solid and dotted segments illustrate the anyon 
matching arising from ii-EC and Zoo-EC respectively. Since 
together they complete a non-trivial loop, the matchings are 
logically inequivalent. 



anyon propagation paths of equal length cross at their 
midpoints, the amplitude corresponding to full propaga- 
tion on both paths can be seen to be at times! It is the 



aim of this section to extend the setting of (Sec. IV B I 



to multiple logical errors while ensuring that the failure 
probability remain exponentially small with system size, 
thereby allowing a fully depolarizing map on the code- 
space with probability exponentially close to 1. 

The basic idea behind this construction is that, for 
large systems, the propagation of the anyons is essentially 
ballistic. Hence, we can divide our lattice into sections, 
and ensure that the paths for anyons of different types 
only cross in regions where we can be (almost) guaranteed 
of the order in which the anyons pass through. It is then 
our task to bound the error probability. 

Let us first consider the probability p s (t) of finding a 
propagated anyon at site s after a propagation time t 
which is given in [42] as 



(t) = f(s; D, sin 2 t) = sin 2s t cos 2 ^^ i , (CI) 



D 



Pi 



where / is the binomial distribution function and D is the 
propagation length (i.e. there are D + 1 possible anyon 
sites in the path). Here time has been normalized such 
that perfect transfer occurs at t = ^. Correspondingly, 
if P is the perfect transfer Hamiltonian for vertex anyons 
and Ho is the projector onto the subspace with a unique 
anyon at the transfer start site, then 



e - ltp n = £ 5 .a s (t)z^n 

where \a s (t) \ = y/p s (t) 



(C2) 



and Z® s is the tensor product of s consecutive Z opera- 
tors along the anyon propagation path. 

We are now in condition to compare an the actual evo- 
lution imposed by two non commuting anyon propaga- 
tions \i/)(t)} an d an ordered idealization of it \ipb;a(t)) 



\xp(t)) = Ue- it( - p °+ p ^W\ip ) 



(C3) 
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We will assume that the physical qubit corresponding to 
the crossing of both paths is between anyon sites s a — 1 
and s a of the anyon path associated to P a and between 
anyon sites s b — 1 and s b of the anyon path associated to 
P b . Furthermore, we will assume s a 3> s b , where what 
is meant by ('3>) will soon be made clear. Under these 
conditions, we will see that \ip(t)) and \ip b ;a(t)) are almost 
equal (at least during the time period corresponding to 
perfect state transfer). 

By definition, we have that (ip(0)\ip b]a (0)) = 1. Let us 
now bound how fast this overlap can actually decay 



d m)\^ b , a {t)) 

d t 



i(ip(t)\[P b ,e 



-itP a 



- itp »tf\i> ). (C4) 



This allows bounding 



d t 



< ||[P b ,e- Jtp °] e " ltp ^|^ )| 



(C5) 



Now let llg' and IL^ be projectors onto the subspace 
with no anyons in the path of P a and the subspace where 
a single anyon is located at the initial site and define II^ 

and analogously. Recalling that |V>o) is a code state 
and our choice of rotation IA, we have 



(nL o) + n< Vl<M = w+l^o) 
(i4 6) + n< 6 Vl<M = u*\fh) 



(C6) 



Commuting these projectors and using the expansion 
( C2 1 of the perfect transfer we may express the RHS of 



equation ( C5 1 by 



5 n, 



(a) iE 



a r (t)X^ r U^U^ ) 



(C7) 



There is only one possible non commuting term in P b 
and this only for s > s a - Furthermore, this term cancels 
for all but two terms in the sum over s' . We may then 



rewrite ( C7 1 as 



2^, J2 s > Sa Vs(t)Z®°Il, 



. . (C8) 

Where J Sb is the strength of the term performing an 
anyon swap between sites s b and s b — 1. Since each coeffi- 
cient accompanies an orthogonal component of the state, 
we may recall the definition in (IC2I) and rewrite ( C8 1 as 



2 j Sb / \p Sb -i(t) + Psb («)] 22 p° (*) sin2 £ 

V s>s a 



(C9) 



where sin 2 e is the amplitude of H^U^U^ \ip<j) . An ex- 
ponentially small upper bound will now be given for the 
expression inside the square root . 

[Ps b -l(t)+Ps b (t)]Es>s a Ps(t) 

< E r <s b nr,D,sm 2 t)Z s > Sa f(s,D,sm 2 t) (CIO) 
= F(s b , D, sin 2 t)F(D - s a ,D, cos 2 t), 



where F(k,N,p) — Yli=of(hN,p) is the cumulative bi- 
nomial distribution function. Assuming 4j < sin 2 t < 4j 



we may use Hoeffding's inequality [43] to bound (CIO I as 



o ' < e D 



(Cll) 



with equality holding for sin t 



2D 



In turn, a tighter 



bound can be obt ained by using Hoeffding's inequality on 
a single factor of (CIO I when sin 2 t > %_or ^ > sin 2 t. 

IV B and s — r > 



1 as in Sec. 



Taking D = NJT 
Sa — Sb ^ D/6 for instance, the obtained upper bound 
becomes exponentially small in N. In turn, the deriva- 



tive ( C5 1 is exponentially small, meaning that the actual 



evolution is approximated by the ordered evolution with 
exponentially good precision in N. 




FIG. 7: In an JV x iV lattice, there are two sets of N/k rows 
(k ~ 0(1)) and two sets of columns and rows, each of which 
corresponds to the construction of (Sec. IV B I for a different 



error type (Xj , Z2 are introduced by columns starting at hor- 
izontal stripes and Z\ and X2 are introduced by rows starting 
from vertical stripes). 

We have formally proven that for two non commut- 
ing anyon propagation paths which intersect with a suf- 
ficiently large offset (i.e. > D/6) the evolution can be 
accurately approximated by ordered anyon propagation. 
There is no obstacle in generalizing this result to many 
such anyon paths, as required to introduce logical errors 
with high probability. Some leading factors of order N 2 
appear but crucial factors remain exponentially decreas- 
ing in N. 

In Fig. [7] we illustrate a configuration allowing the 
simultaneous introduction of all possible logical errors 
by anyon propagation within the lattice. The marked 
stripes of width N / k indicate locations where perpendic- 
ular anyon propagations begin or end. The perturba- 
tion to be introduced is chosen randomly as in Sec. |IV B| 
such that each propagation row/column starts with equal 
probability in either of each pair of opposing stripes. Tak- 
ing k fixed allows sufficiently many repetitions of the sin- 
gle row/column construction that the probability of in- 
troducing each type of logical error approaches g expo- 
nentially fast with N. Thus, after a perturbed evolution 
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for time tf = 0(N), and a final application of an ar- tially close to the maximally mixed state 1/4 of the code- 
bitrary error correcting protocol based on unperturbed space, 
syndrome measurement, the resulting state is exponcn- 



[1] P. W. Shor, SIAM J. Comput. 26, 1484 (1997). 
[2] A. M. Childs and W. van Dam, e-print arXiv:0812.0380 
(2008). 

[3] A. W. Harrow, A. Hassidim, and S. Lloyd, Phys. Rev. 

Lett. 103, 150502 (2009). 
[4] C. H. Bonnet and G. Brassard, in Proceedings of the 
IEEE International Conference on Computers, Systems 
& Signal Processing, Bangalore, India (IEEE, New York, 
1984), pp. 175-179. 
[5] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 
[6] P. W. Shor, Phys. Rev. A 52, R2493 (1995). 
[7] D. Aharonov and M. Ben-Or, quant-ph/9611025 (1996). 
[8] D. Gottesman, Phys. Rev. A 57, 127 (1998). 
[9] D. Gottesman, e-print arXiv:0904.2557 (2009). 
[10] A. Y. Kitaev, Annals of Physics 303, 2 (2003). 
[11] E. Dennis, A. Kitaev, A. Landahl, and J. Preskill, Jour. 

Math. Phys. 43, 4452 (2002). 
[12] T. Kato, Perturbation Theory for Linear Operators 

(Springer- Verlag, Berlin, 1966). 
[13] M. B. Hastings and X. G. Wen, Phys.Rev. B 72, 045141 
(2005). 

[14] A. Hamma and D. A. Lidar, Phys. Rev. Lett. 100, 030502 
(2008). 

[15] R. Alicki, M. Fannes, and M. Horodecki, e-print 

arXiv:0810.4584 (2008). 
[16] R. Alicki, M. Horodecki, P. Horodecki, and R. Horodecki, 

e-print arXiv:0811.0033 (2008). 
[17] S. Chesi, D. Loss, S. Bravyi, and B. M. Terhal, e-print 

arXiv:0907.2807 (2009). 
[18] F. Pastawski, A. Kay, N. Schuch, and I. Cirac, Phys. Rev. 

Lett. 103, 080501 (2009). 
[19] S. Bravyi, D. P. DiVincenzo, D. Loss, and B. M. Terhal, 

Phys. Rev. Lett. 101, 070503 (2008). 
[20] M. M. Wolf, Nat. Phys. 4, 834 (2008), ISSN 1745-2473. 
[21] D. Bacon, Phys. Rev. A 78, 042324 (2008). 
[22] S. Chesi, B. Rothlisberger, and D. Loss, e-print 



arXiv:0908.4264 (2009). 
[23] A. Hamma, C. Castelnovo, and C. Chamon, Phys. Rev. 

B 79, 245122 (2009). 
[24] D. Bacon, Phys. Rev. A 73, 012340 (2006). 
[25] A. Kay, Phys. Rev. Lett. 102, 070503 (2009). 
[26] C. King, Trans. Inform. Theory 49, 221 (2003), ISSN 

0018-9448. 

[27] P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 
[28] L. C. Venuti and P. Zanardi, Phys. Rev. Lett. 99, 095701 
(2007). 

[29] D. Kribs, R. Laflamme, and D. Poulin, Phys. Rev. Lett. 

94, 180501 (2005). 
[30] M. Christandl, N. Datta, A. Ekert, and A. J. Landahl, 

Phys. Rev. Lett. 92, 187902 (2004). 
[31] P. Karbach, Phys. Rev. A 72, 030301 (2005). 
[32] A. Kay, e-print arXiv:0903.4274 (2009). 
[33] S. Bravyi and B. Terhal, New J. Phys. 11, 043029 (2009). 
[34] A. Kay and R. Colbeck, e-print arXiv:0810.3557 (2008). 
[35] P. R. A. Vaidya, Geometry helps in matching (ACM, 

1988), pp. 422-425, ISBN 0-89791-264-0. 
[36] H. L. Richards, S. W. Sides, M. A. Novotny, and P. A. 

Rikvold, Journal of Magnetism and Magnetic Materials 

150, 37 (1995), ISSN 0304-8853. 
[37] E. Cirillo and J. Lebowitz, Journal of Statistical Physics 

90, 211 (1998). 
[38] G. Grinstein, IBM J. Res. & Dev. 48 (2004). 
[39] K. R. Brown, Phys. Rev. A 76, 022327 (2007). 
[40] D. A. Lidar, Phys. Rev. Lett. 100, 160506 (2008). 
[41] S. P. Jordan, E. Farhi, and P. W. Shor, Phys. Rev. A 74, 

052322 (2006). 

[42] A. Kay and M. Ericsson, New Journal of Physics 7, 143 
(2005), ISSN 1367-2630. 

[43] W. Hoeffding, Journal of the American Statistical Asso- 
ciation 58, 13 (1963), ISSN 01621459. 



